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ABSTRACT 

The canopy-air temperature difference (6T) may provide an index for scheduling irrigation. 
Combining the Monteith transpiration equation with both uptake from a single-layered root zone 
and change in internal storage of the plant, we have explicitly solved the continuity equation for 
water flux in the soil-plant-atmosphere system. Using appropriate parameters for com (Zea mays L.) 
the model indicates that both daily total transpiration and soil-induced depression of plant water 
potential may be inferred from mid-day 6T. For the soil-plant-weather data used in the simula- 
tion, either a mid-day spatial variability of about 0.8K in canopy temperatures or a field-averaged 
6T of 2 to 4K might be a suitable criterion for irrigation schedul ig. 
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INTRODUCTION 

Irritation scheduling should maximize the benefit of each unit of irrigation water. Recently 
Jackson (1982) argued that direct observation of some plant factor should be a superior approach to 
optimizing irrigation than use of soil or meteorologic factors. The plant factors that have been 
studied as indicators of plant stress are leaf water potential, stomatal resistance, leaf temperature, 
and canopy temperature. The leaf-based measurements (potential, resistance, and leaf temperature) 
are time consuming when fields of the size of hectare or larger are to be sampled to obtain a repre- 
sentative value. Therefore, according to Jackson (1982), the most promising approach to irrigation 
scheduling would be measurements of canopy temperature using infrared radiometers. Further- 
more, identification of crop stress using infrared radiometers can be done at large scale using space- 
borne sensors. 

Canopy temperature represents an integrated response of a crop to prevailing weather and 
soil water conditions. In the absence of any soil water deficit, diurnal and day-to-day variations 
in the canopy temperature would be due to the weather parameters (insolation, wind speed, air 
and dew point temperatures). As transpiration becomes limited by depleting soil water the canopy 
temperature increases with respect to air temperature, so that the reduced transpiration is balanced 
by increased heat loss. In a field study of irrigation scheduling for com ( Zea mays L.) using infra- 
red thermometers, Qawson and Blad (1982) showed that both canopy temperature and its spatial 
variation can be used effectively to reduce irrigation water without affecting yield. 

The plant, soil, and atmospheric factors that control plant water stress and leaf-air temperature 
difference are fairly well-known and can be mathematically expressed. Therefore, simulation pro- 
vides a way to examine crop behavior without the labor and uncontrolled aspects of field experi- 
ments. 

In this paper we develop a simulation model for canopy-air temperature difference for a com 
crop and use it to study temporal and spatial variation of the temperature difference with respect 
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to irrigation. Choudhury (1983 a,b) solved Monteith’s (1965) and van den Honert’s (1948) equations 
for transpiration and tested the results against Idso’s (1982) data for unstressed com and soybeans 
(Glycine max L.) on clear days. The continuity equation for water flux in the soil-plant-atmosphere 
system is solved by expressing the stomatal resistance and the rate of plant tissue water loss in terms 
of a plant water potential (Federer, 1 979). We have also added the dynamics of root-zone soil water 
following Feddes and Rijtema (1972), and the increase in flow resistance in a drying plant (Boyer, 
1971 ; Nulsen and Thurtell, 1978). The model currently ignores soil • vaporation. Dynamic equations 
then show the progression of soil water and plant water storage and their effect on canopy temper- 
ature through the energy balance equation. 


By examining the temporal variation of canopy-air temperature difference during drying we 
determine if it can be used for irrigation scheduling. We also simulate expected variation of canopy 
temperature across a field to determine its effect on scheduling. 

DESCRIPTION OF THE MODEL 
Vapor Phase 

For transpirational flux, E A , (m/s), we use Monteith’s (1981) modification of Monteith’s 
(1965) combination equation. This modification accounts for the dependence of net radiation on 
the canopy temperature by defining an effective boundary-layer resistance. 

C — ^ * ^~P Pm ( e * ~ e « )/ r t /n 

A Ly (A+ 7 (r, +r c )/r t ] 

where R no is the net radiation flux absorbed by the canopy if the canopy surface were at air tem- 
perature (W m -3 ), e* and e # are, respectively, the saturated vapor pressure (kPa) at air and dew point 
temperature, A is the slope of saturated vapor pressure at the air temperature (kPa K -1 ), C p and 
p t are, respectively, the heat capacity (J kg' 1 K" 1 ) and density (kg m* 3 ) of air, y is the psychro- 
metric constant (kPa K" 1 ), L r is the latent heat of vaporization (J m -3 ), r c is the canopy resistance 
(s nf 1 ), and r ( is the effective boundary -layer resistance for heat and longwave radiative transfer 
(s m" 1 ), given by Monteith ( 1 98 1 ) as 
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( 2 ) 


r *f\ + 

where e c is the canopy emissMty, a is the Stefan-Boltzman constant, T a is the air temperature and 
r a is the usual boundary -layer resistance (s m~ l ), calculated according to Grace (1977) for com 

75 

r « * LAI U 1 '* 

where LAI is the leaf area index and U is the wind speed (m s’ 1 ). The net radiation at the top of 
the canopy is calculated from global insolation (S, W m' 3 ), albedo (a), longwave emissivities of 
canopy (e c ) and air (e a ) (Idso, 1 98 1 ) and air temperature as 

R B| -(l-o)S + (e a -e e )oT^ (4) 

From the observed (Impens and Lemeur, 1 969) extinction coefficient of net radiation in a com 
canopy, the absorbed net radiation is 

R „o •R.ill- exp { 0.055 LAI 3 - 0.622 LAI | ] (5) 

The canopy resistance depends most strongly upon leaf-water potential and global insolation 
(Kramer, 1969). From the observed insolation dependence by Uchijima (1976) and leaf-water 
potential dependence by Reicosky and Lambert (1978) we have calculated the canopy resistance 
from 


( 4T$7 ) 


1.4 X 10 4 ft 175, 

r * [ 1 + 1 

c LAI(*-* e ) S+10 

where \p is the plant water potential (m) and \p c is the critical potential for stomatal closure. 


( 6 ) 


Liquid Phase 

The rate of change of plant water storage (E T ) is modeled according to Federer (1979). 
Milne et al. (1983) recently justified this approach for trees. The stored water, Q (mm), in the 
plant is assumed to be at a potential \p^ that differs from the plant water potential, ip, which 
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controls stomatal closure and water flux through the plant. The stored water is considered to be 
isolated from the transpiration stream by a resistance R q (s). Thus, 


Et 



( 7 ) 


Federer (1979) modeled as a linear function of Q 


( 8 ) 


where Q o is the maximum storage. However, for the current model we tried to get a more realistic 
function. Presumably the stored water is in living cells, primarily in the plant stem; but we are not 
aware of - Q relations ot such cells. As the best alternative we assume that the relation of leaf 
water content to leaf water potential applies also to storage cells in the stem. To observations of 
Reicosky and Lambert (1978) for com we fitted a bilinear equation (Fig. 1) 


where 


- -450 Xf - ( 1 3 1 X + 84)(l-f) 



(9) 


and 


1 

1 + 6093 X 5 


The rate of water extraction by plant roots (E $ in m s -1 ) is calculated from van den Honert’s 
1 948 equation 


' R . +R p 


( 10 ) 


where is the soil water potential (m), and R ( and R p are, respectively, the resistances for water 


flow within soil to the root surface and from the root surface to leaf stomata. 
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The soil resistance is calculated from an empirical form of the Gardner - Cowan equation 
given by Feddes and Rijtema (1972) 


R. 


0.0013 

" Z.fr *<*.> 


(ID 


where Z #ff is the effective root depth (m) and K(^ t ) is the unsaturated hydraulic conductivity of 
soil (m s'* ), given in the parametric form by Campbell (1974) as 

J. S ♦ J/b 

K(* ( )«Kto (12) 

where the parameters K , 4* mt and b for different soil texture classes can be found in Gapp and 
Homberger (1978). 


The plant resistance is calculated as the ratio of root resistance per unit length r, (s m _1 ) and 
root length per unit area L,. (m“* ) as 


R 


p 



(13) 


Boyer (1971) and Nulsen and Thurtell (1978) observed that when a com plant is subjected to stress 
its root resistance increases. Because shrinkage of roots resulting from a loss of tissue water is 
thought to increase the resistance (Vaadia et al., 1961), we have modeled the stress induced changes 
in the plant resistance as 

r r - r^ { 1 1 - f 1 X” + l| (14) 

where r ra is the unstressed root resistance, and f and X are tissue water dependent quantities de- 
fined in Equation (9). The functional form is chosen quite arbitrarily; however, Equation (14) 
recognizes a non-linear dependence of the resistance on the storage potential, , and on the dura- 
tion of stress noted by Nulsen and Thurtell (1978), as reflected in the amount of storage. 


Equarioru of Continuity and Energy Balance 

The continuity equation for water flux in the soil-plant-atmosphere system is 


5 




ORIGINM- PAGE » 

„ roo* QUAEIW 

E a -E, + E t (15) 

where the fluxes are defined by Equations (1), (7) and (10). 

The energy balance equation is 

T.-T.-tR.-L v E A ](r -f /C p p.) (16) 

where T c is the canopy temperature. 


Dynamic Equations 

The dynamics of plant storage water Q (mm) follows from Federer (1979) 

7^- -10 s E t 
dt T 


(17) 


and the dynamics of root-zone soil water potential comes from integrating the diffusion equa- 
tion for moisture flux 




(18) 


where the soil capacitance V(^) is given by (using Clapp and Homberger’s (1978) equation relating 
volumetric moisture to the matric potential) 


d 6 0 . i*l /b 

V(* t )«(— ) — r-f- (f 51 -) 

• d^, b + mt \p t 


(19) 


and 6 % and are, respectively, the saturation volumetric moisture and matric potential. 


Method of Solution 

Substituting ( 1 ), (6), (7), and (10) into (15) yields a quadratic equation in leaf water poten- 
tial which has the solution 

4> - [- B + (B 2 - 4 AC)* ] /(2A) (20) 

where 
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B • (A + 7*) r, (E 0 -*,)-(«- it') tj-i- + ^-l 


C * (a - 7 r ) \l/ e - E 0 a 

7 * m 7Tjr t 

c .^R»o + C p p.(c;-e.)/r t 

E - — wu + 7*i 

* - *• t ♦« 

• R. + R, R, 

a * (A + 7 *) r, \lf e 


1 4 X 10 4 

r' — [1 + 175/(8+10)1 


For any given set of weather variables (insolation, wind speed, and air and dew point temperatures), 
plant properties (leaf area index, storage resistance, root depth and root length per unit area) and 
soil texture parameters (b, K , ^ g and 0 ( ) one can calculate the plant potential when the plant 
storage Q and the soil water potential are specified. Diurnal weather variables determine the 
evaporative demand (E o ) and Equations ( 1 7) and ( 1 8) give the diurnal variation of \p t and Q. 

The numerical method used for the solution of Equations ( 1 7) and (18) is the Douglas-Jones 
implicit linearization (Douglas and Jones, 1963). If we define functions F, (0,, Q, \l>) and F 2 


(0,.Q. tfO* s 




F, <*,,<}, *)— 10* &!*] 
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then from the starting values and Q* one obtains the new values and Q" after time At in 
two steps; the intermediate step, 

J,-*; ♦yp. (23> 

3-<r (24) 

and the final step 

+ At F, (^,5, f) (25) 

Q" "Q 0 +AtF a (J t ,5.?) (26) 

where \l/° and are obtained from Equation (19) using, respectively, (\if° , Q° ) and (^ s , $). 

Knowing ^ from Equation (20) one can calculate the transpirational flux E A from Equations 
(6) and (1). The canopy temperature T c then follows from Equation (16). 

METHODS 
Panmeter Selection 

Reicosky et al. (1975) and Reicosky and Lambert (1978) provide concurrent Held data for 
microclimate and leaf water potential for corn crops. We chose data from June 1 and June 6 in 
Reicosky et al. (1975) and Day 167 in Reicosky and Lambert (1978) for testing and running our 
model. Effective rooting depth (Z #ff ) for these crops was 0.3 m; the leaf area index (LAI) was 3.5 
on June 1 , 3.6 on June 6, and 4.2 on Day 167. We assume albedo (a) is 0.22, canopy emissivity 
(< c ) is 0.97, and critical potential for stomatal Mosure (^ c ) is -2 1 5 m. For unstressed root resist* 
ance (r ro ) we take Newman’s (1973) value for detopped com. 4 X 10 ,J s m*‘ . Fro-:'. wet and 
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dry biomass data of com (Want ct al. f 1 982) we estimate the maximum plant water storage (Q # ) 
is 2 mm. 


From Clapp and Homberger (1978) we chose soil parameters for a sandy loam soil: b ■ 4.9, 

6 % - 0.435, - - 0.22 m, K, ■ 3.5 X 10"* m s“* . 

We chose root length per unit area (L T ) and storage resistance (R q ) as the values that best fit 
the model leaf water potential to the field data (assuming midnight soil water potential of - 0.6 m): 
(Lj. , R q ) values are (5 X 10 s , 4 X 10 l# ) for June 1,(4 X 10 s , 4 X 10 l# ) for June 6, and (5 X10 s , 

4 X 10*) for Day 167. Good agreements shown in Fig. 2 provide some confidence in the simula- 
tion model, although they do not provide an ab initio verification. A sensitivity study, in which 
the microclimatic data are varied without changing the crop parameters or varying the crop param- 
eters without changing the microclimate, is not done in this paper. 

Field-Averaged Drying Cycle and Irrigation 

Starting with a mid-night soil water potential of -0.6 m we simulated soil drying for ten days 
with soil, crop, and climate data of June 1 , June 6, and Day 167. This ten-day drying period is 
refened to below as the drying cycle. Following this drying cycle a simulated irrigation was exe- 
cuted by setting the soil water potential to -0.6 m on the 0100 hour of the eleventh day of simula- 
tion. A second ten-day drying period was then simulated. 


Spatial Variation 

The spatial variation of soil hydraulic properties is assumed to be described by scaling theory 
(Philip, 1981). According to this theory stochastic variation of saturated conductivity and potential 
are described by the equations 
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(27) 

( 28 ) 


k,.«»K, 

* V«/t 

where K and J m{ are the field averafed values and the scaling parameter t is a random number. 

We assumed t is a log-normally distributed random number with mean -0.139 and standard devia- 
tion 0.5 1 1 , after Warrick et al. ( 1 977) for Panoche soil. Ten random scaling parameter values were 
generated, and for each of the ten sets of hydraulic parameters and each of the three test days we 
simulated a ten-day drying cycle. 

We also considered spatial variation of applied irrigation water. After each of the above drying 
cycles the amount of applied water was generated stochastically by assuming a uniform distribution 
with mean equal to the amount of water depleted from the root-zone and a coefficient of variation 
of 0.1. The choice of the coefficient of variation is according to the observations of Clawson and 
Blad (1982) for a sprinkler irrigated field. Following irrigation, at 0100 hour, we simulated a 
second ten-day drying cycle. This sequence of ten drying cycles consider spatial variability of 
applied irrigation water in addition to the soil hydraulic heterogeniety. 

SIMULATION RESULTS AND DISCUSSION 
Field-Averaged Results 

For all three data sets the canopy temperature at 1300 hour increases fairly slowly at the 
initial stages of soil drying, then the later increase fairly rapidly (Fig. 3). This trend in canopy tem- 
perature is consistent with the observations of Clawson and Blad (1982) who had noted that in 
sandy soils once stress develops, it progresses fairly rapidly. The diurnal trends of canopy-air tem- 
perature differences are shown in Fig. 4. 

The canopy temperatures do not return to their original unstressed values after irrigation 
(Fig. 3). This is consistent with the observation: of Jackson et al. (1981). The canopy tempera- 
tures continue to decrease for a few days after irrigation before beginning to increase again as the 
soil continues to dry. This trend in the simulation appears because with irrigation, the incurred 
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km of plant tissue water during the drying period begins to get restored by the soil. Due to con- 
tinual soil drying, however, a complete restoration of the plant tissue water does not occur; and, 
therefore, the canopy temperatures never return to their original (Day 1) unstressed values. It was 
found that when irrigation is done after 5 or 6 days of drying, then the canopy temperatures do 
return to values close to the day 1 of simulation. The leaf water potentials and stomatal resist- 
ances also show a recovery period of a few days after irrigation (Fig. 5). Such a lag in the recovery 
of ieaf water potentials has been observed for several crops (Kramer, 1969). The rate of recovery 
of the leaf water potential is more rapid on high evaporative demand day (Day, 167) than on low 
evaporative demand days (June 1). 

The recovery lag suggests that if the stress history of a crop is not known, it may not be 
possible to infer unambiguously irrigation needs from measuring canopy temperature, leaf water 
potential or stomatal resistances. For example, Figure 3a shows that day 1 and day 1 1 have almost 
identical “well-watered” soil water potentials, but the canopy-air temperature difference differs 
by about 0.8K. Plant storage, Q, recovers slowly both because of its own resistance R q , and because 
the root resistance r r depends on Q. Consequently the plant may take one to several days to 
respond to changes in soil water. 

Clark and Filer (1973) compared leaf water potential, stomatal resistance and leaf-air tempera- 
ture difference of well-watered and water-stressed peas in order to determine the most suitable crop 
water stress indicator. They concluded that the leaf water potential was more responsive to changes 
in plant water status. The relationship of the leaf water potential to the soil water potential and 
the canopy-air temperature difference are shown in Figures 6 and 7. The relationship between the 
two potentials is essentially linear, for any particular day, down to a soil water potential of about 
-60 m, below which the plant potential decreases less rapidly (Fig. 6). A linear relationship 
between the two poter -* 1 ; has been observed by Shinn and Lemon (1968), Ehling et al. (1968) and 
Brady et al. (1974), among others. From Figure 7 one would conclude that if leaf water potential 
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is • good indicator of crop water status, so should be the canopy-air temperature difference. How- 
ever, when simulation results for all 3 days are composed to form one “data” set, it is difficult to 
see how either the leaf water potential cr the canopy-air temperature difference could be a reliable 
indicator of irrigation needs. In Fig. 5 a leaf water potential of-175 m would correspond to a soil 
water potential of about -65 m for June 1 , but the same leaf water potential for Day 167 corre- 
sponds to a soil water potential of about —20 m. Similarly, from Fi®s. 3a and 3b we see that a 
canopy-air temperature difference of 3.5K corresponds to a soil water potential of -45 m for June 
1 , and -70 m for Day 167. 



The atmospheric variability significantly affects both the leaf water potential ($) and the 
canopy-air temperature difference 5T (Idso et al., 1981 ; Jackson et al., 1981). At high soil water 
potentials if we neglect the plant storage change and the soil resistance for water flow to the root 
surface, then equations (10), (15), and (16) lead to 

I r 

*T-[R„+=^-] (—M (29) 

Rf C p P . 

where the superscript o is used to indicate high soil water potentials (the day 1 values are assumed). 
Idso (1982) has shown that 5T° is a linear function of atmospheric vapor pressure deficit, and 
that it is relatively insensitive to all other weather parameters. Then equation (29) shows that i//° 
will also be strongly affected by the atmospheric vapor pressure deficit, through a linear relation- 
ship. Observations of Idso et al. ( 1 98 1 ) confirm such a linear relationship. Jackson et al. ( 1 98 1 ) 
defmed a crop water stress index (CWSI), 


CWSI- 


5T- 8T° 



( 30 ) 


to normalize weather dependent variability of canopy temperature. We can define an analogous 


normalization for the plant water potential, a plant water stress index (PWSI), as 
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(3l) 

where is the critical potential for stomatal closure. Figure 8 shows that PWS1 normalizes for the 
atmospheric variability and that the soil water information can be inferred more-or-less unam- 
biguously via this index. The relationship between PWSI and CWSI is shown in Figure 9. If radio- 
metric remote sensing is used to calculate a CWSI (Jackson et al., 1981 ) then one can also obtain a 
PWSI and hence, the soil water potential. 

The concept of using aircraft and space-bome radiometric data to predict crop water require- 
ments for large areas has been investigated by Bartholic et al. (1972), Schmer and Werner (1974) 
and Jackson et al. (1977), among others. One aspect of this investigation has been to calculate 
large area daily evapotranspiration (ET) from one-time-of-day measurement of ground temperature. 
Jackson et al. (1977) observed a significant linear relationship between daily ET and post noon 
canopy-air temperature difference for wheat. The present simulation shows that for individual days 
the daily total transpiration is almos. linearly related to the 1300 hour canopy-air temperature 
difference, and the results for the 3 days taken together may also be approximated by a linear rela- 
tionship (Figure 10). Comparison of these simulation results with the observations of Jackson 
et al. (1977) indicates that these linear relationships would at least be crop-specific, since for wheat 
a canopy-air temperature difference of 2K corresponded to an ET of about 3 mm, whereas we 
simulate 6 mm for com. The validity of the present simulation results remains to be confirmed by 
observations. 

Spatial Variability 

The range of canopy temperatures (CTV) at the 1 300 hour of each day of the drying period 
are shown in Figures 1 la and b. The results in Figure 1 la are from the first drying cycle which 
considered soil hydraulic heterogenie ty, and in Figure 1 lb from the second drying cycle which 
started following the above drying cycle by stochastically replenishing the amount of depleted soil 
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water from the root-zone. The simulated steep increase in CTV after the first several days of drying 
is fully consistent with the observations of Clawson and Blad (1982) and supports their contention 
that CTV can be used very effectively for triggering irrigation. In Clawson and Blad’s field experi- 
ment, irrigation was triggered at CTV ■ 0.8, and we see in Figure 1 la and b that this occurs first 
for Day 167, then June 1 and finally June 6. This arrangement of days follows the daily total 
transpiration (Figure 10), and hence the rate of soil drying. The onset of CTV ■ 0.8 occurs about 
a day earlier in the second drying cycle (Figure 1 1 b) as compared to the first cycle (Figure 10a). 

An early onset of an equivalent stress condition in the second drying cycle would have been ex- 
pected because this drying cycle got started with a stressed crop. A spatial variability of irrigation 
water in addition to the soil hydraulic heterogeniety leads to a little larger spatial variability of 
canopy temperatures (Fig. 1 1). 

In order to quantify the soil-water induced stress associated with the spatial variation of 
canopy temperatures, we plotted the CTVs from Figure 1 la as a function ni P9/SI (Figure 12). A 
CTV * 0.8 corresponds to PWSI between 0.3 and 035, which, from Figure 8, give the root-zone 
soil water potential between -30 and -35 m or a depletion of conventionally defined available 
water between 65% and 70%. For the soil-plant-atmosphere data used in the present simulation, we 
are essentially substantiating the long-held wisdom of irrigation scheduling, albeit via the modem 
technology of infrared remote sensing. (For the CTVs in Figure 1 lb, however, the crop stress con- 
ditions cannot be related unambiguously to the soil water status because of the lag in recovery of 
canopy temperatures.) 

For irrigated com crops Sumayao et al. (1982) found the leaf stoma tal resistances to be less 
than 500 s m* 1 , which, according to the present simulation, is consistent with CTV - 0.8 or soil- 
water potentials between -30 and -35 m (cf., Figures 3 and 5). The field-averaged canopy-air 
temperature differences are between 2 and 4K when the soil water potentials are between -30 and 
-35 m. 
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Our simulations indicate that either a 2-4* canopy-air temperature difference or about 1*C 
spatial range in canopy temperatures at midday might be suitable indicators of need for irrigation 
in com. 
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CAPTION TO THE FIGURES 

Figure 1 . Relationship of plant storage water potential and plant water content. See text for the 
meaning of symbols. 

Figure 2. Observed microclimatic data, and observed and simulated leaf water potentials. The 
data in (2a) and (2c) are from Reicosky et al. (1975), and (2b) from Reicosky and 
Lambert (1978). 

Figure 3. Simulated canopy-air temperature difference and soil water potential at 1 300 hour of 
different days. A simulated irrigation is done at 0100 hour of the 11th day. The air 
temperature is constant for each data set, so variation in the canopy-air temperature 
difference is caused by variation in canopy temperature. 

Figure 4. Dirunal variation of canopy-air temperature difference. The curves are labeled by the day 
of simulation as shown in Fig. 3. 

Figure 5. Leaf water potential and leaf stomatal resistance at 1300 hour of various days of simula- 
tion. Note the lag in recovery after irrigation at 01 00 hour on the eleventh day. 

Figure 6. Relationships of 1300 hour leaf and soil water potentials for the three data sets. 

Figure 7. Relationships of 1300 hour leaf water potential and canopy-air temperature difference 
for the three data sets. 

Figure 8. Relationship of 1 300 hour plant water stress index as defined by equation (3 1 ) and soil 
water potential. 

% 

Figure 9. Relationships of 1 300 hour plant water stress index and crop water stress index. 
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Figure 10. Relationships of daily total transpiration and 1300 hour canopy -air temperature 
difference. 

Figure 1 1 . The range of canopy temperatures at 1 300 hour for the three data sets, (a) the 

variability originating from spatial variation of soil hydraulic properties, (b) the vari- 
ability originating from spatial variation of soil hydraulic properties and irrigation 
water application. 

Figure 1 2. Relationships of 1300 hour plant water stress index and spatial variability of canopy 
temperatures resulting from soil hydraulic heterogeniety. 
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